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NOMENCLATURE 
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D 
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®w 

n 


<lw 

Ow 

v a 

v w 

v d 


Frontal area of nacelle, (m^) 

Droplet Diameter (mm) 

Liquid Water Content (gm/m 3 ) 

Mass flow rate of water (gm/sec) 

Number of drops counted 

Total liquid volume in the measurement volume (cm 3 ) 
Volumetric water ingestion rate (cm 3 /sec) 

Airspeed (m/sec) 

Water jet speed at the spray nozzle (m/sec) 

Average water drop speed (m/sec) 

distance downstream from the spray nozzle (m) 


Greek Letters 

p w Water density (gm/cm 3 ); 1.0 gm/cm 3 
5 Thickness of the Light Sheet (mm) 


Subscripts 

a Average spray property 
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EXECUTIVE SUMMARY 


Water droplet ingestion into turbine engines resulting from heavy rain and 
wheel spray generated on a wet runway is of importance. Adverse effects of 
large quantities of water ingestion can include the compressor stall and 
combustor flame-out. Engine certification requirements as set forth in FAA 
regulations call for continued engine operation at takeoff and flight idle 
conditions while ingesting water at 4 percent by weight of airflow, generated 
by a spray to simulate rain. There is also a certification requirement on the 
entire aircraft system which dictates that the system must be designed to 
prevent hazardous quantities of water from being ingested into the engine 
during takeoff, landing and taxiing operations on wet runways. The present 
work was undertaken to develop measurement techniques of two -phase droplet 
laden airstreams during engine water ingestion. The ultimate objective is to 
correlate the non- intrusive measurements of the water ingestion rate and 
droplet size and spatial distribution at the engine inlet with engine 
performance parameters. Such techniques and data will assist the FAA in 
evaluating current water ingestion certification tests. 

A non- intrusive optical technique was developed for the determination of liquid 
mass flux in a droplet laden airstream. The technique is also capable of 
providing information on the droplet size and spatial distribution at the 
nacelle inlet plane. 

Independent measurements of the liquid water content (LWC) of the droplet laden 
airstream and of the droplet velocities were made at the inlet plane of a 
simulated nacelle in a wind tunnel for the liquid mass flux determination. The 
liquid water content was determined by illuminating and photographing the 
droplets contained within a thin slice of the flow- field by means of a sheet of 
light from a pulsed laser. Fluorescent dye introduced in the water enhanced 
the droplet image definition. The droplet velocities were determined from 
double exposed photographs of the moving droplet field. The technique was 
initially applied to a steady spray generated in a wind tunnel. It was found 
that although the spray was initially steady, the aerodynamic breakup process 
was inherently unsteady. This resulted in a wide variation of the instantane- 
ous liquid water content of the droplet laden airstream. The standard 
deviation of ten separate LWC measurements was 31 percent of the average. 
However, the liquid mass flux calculated from the average LWC and droplet 
velocities came within 10 percent of the known water ingestion rate. 
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INTRODUCTION 


The effects of water droplet ingestion into turbine engines resulting from 
heavy rain and wheel spray generated on a wet runway is of importance. The 
effects of water ingestion on engine performance have recently been 
investigated (reference 1) and a probe for stagnation pressure measurement in a 
droplet laden airflow was developed (reference 2) . The adverse effects of 
large quantities of water ingestion can include the compressor stall and 
combustor flame-out (reference 1). Engine certification requirement as set 
forth in FAA regulations (reference 3) call for continued engine operation at 
takeoff and flight idle conditions while ingesting water at 4 percent by weight 
of airflow, generated by a spray to simulate rain. There is also a 
certification requirement on the entire aircraft system which dictates that the 
system must be designed to prevent hazardous quantities of water from being 
ingested into the engine during takeoff, landing, and taxiing operations, 
(reference 4) . 

The criteria for both water ingestion certification tests are somewhat 
arbitrary. The 4 percent by weight water ingestion test for engine 
certification does not address such issues as droplet size and their spatial 
distribution over the frontal area of the nacelle. For performance on a wet 
runway, the criterion is even more arbitrary in that it does not specify what 
constitutes a hazardous quantity of water ingestion. 

The present work was undertaken to develop measurement techniques in two-phase 
droplet laden airstreams to better quantify the engine water ingestion. The 
ultimate objective is to correlate non- intrusive measurements of the water 
ingestion rate and droplet size and spatial distribution at the engine inlet 
with engine performance parameters. Such techniques and data will assist the 
FAA in evaluating current water ingestion certification tests. 


REVIEW OF APPLICABLE MEASUREMENT TECHNIQUES 

A review of available measurement techniques for sprays was undertaken to 
determine the applicability of the existing techniques to quantitative 
determination of water ingestion into an engine. The quantities that need to 
be determined are: 

a) Drop size distribution together with spatial distribution of drops at 
the nacelle inlet. 

b) The mass flow rate of water crossing the nacelle inlet plane at any 
instant. 

Extensive work on measurements of wheel sprays generated by aircraft 
under -carriages was carried out by Barrett (reference 5). A spray intensity 
probe was developed to measure mean local dynamic pressure generated by moving 
droplets. However, all of Barrett's measurements were in the near field of the 
wheel generating the spray and the technique used only provided time averaged 
local measurements in the spray. 

Several mechanical, electrical, and optical methods are available for droplet 
size determination in fuel sprays as surveyed in review articles by Jones 
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(reference 6) and, McCreath and Beer (reference 7) . None of the techniques 

surveyed would be capable, in their existing form, of satisfying the second 
requirement above. It was considered possible that the requirement could be 
met by some modification of the existing techniques . 

The first trial approach was an extension of the charged wire probe technique 
described by Gardiner (reference 8) for drop size determination in water 
sprays. In that technique, the electrical pulse generated in a circuit 
containing a charged electrode when a drop impacts the electrode is measured by 
a pulse height analyzer. The probe is initially calibrated using known size 
drops impacting the electrode. A relationship is established between the pulse 
height and the drop size. In practice, the drop size characteristics are 
derived from the pulse height statistics stored in a pulse height analyzer. 
The technique as described by Gardiner (reference 8) thus provides the local 
drop size distribution in a spray. In the present work, a modification of the 
technique was considered. Instead of a single electrode probe, a charged grid 
was considered. A copper wire mesh with wire spacing of approximately 1 mm was 
charged to 2000 volts by a high voltage DC power supply. The idea was to 
install the screen at the nacelle inlet plane so that all drops larger than the 
mesh size would be intercepted by the screen. The charge transfer between the 
screen and the impacting drops would set up a current in the circuit supplying 
the grid. The current would be proportional to the total surface area of the 
drops impacting per unit time. Then if the size distribution were to be 
determined by an independent optical technique, the volume flow rate of water 
would be proportional to the product of the Sauter Mean Diameter (SMD) of the 
spray and the current supplied to the grid. 

In practice, however, several difficulties were encountered with this 
technique . 

1) The charge distribution on the screen was non-uniform, resulting in 
different pulse characteristics for different impact locations on the 
screen for the same size drop. 

2) Wetting of the screen altered the initial charge distribution. 

3) Water film on the support set up a conduction path from the screen to 
the ground, causing a leakage current. 

Because of these difficulties and the fact that an optical technique was still 
needed in conjunction with the charged screen for water flow determination, 
this approach was abandoned in favor of a purely optical non- intrusive 
technique . 


PRESENT TECHNIQUE 

A non- intrusive optical technique was developed for the determination of drop 
sizes, spatial distribution of drops at the nacelle inlet plane and 
instantaneous mass flow rate of liquid water entering the nacelle. The liquid 
water mass flow rate is determined by independent measurements of the liquid 
water content of the droplet laden airstream and the droplet velocity at the 
nacelle inlet plane. 
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For liquid water content (LWC) determination, a thin cross-section of the flow 
close to the nacelle inlet plane is illuminated by a pulsed laser light sheet. 
The drops contained within this light sheet are photographed by a camera placed 
with its axis nearly normal to the plane of the light sheet (figure 1) . The 
duration of illumination of drops by the laser pulse is extremely short, about 
10 nano second. Therefore, the motion of the drops is frozen in the 
photographs . The thickness of the laser sheet is controlled by a beam expander 
and a slit. The slit was of 9 mm width while the initial beam width was 14 mm. 
The slit thus allowed only the intense central portion of the expanded beam to 
pass through, cutting off the less intense outer portions. 

The depth of field of the camera is set to be larger than the light sheet 
thickness by proper selection of aperture and magnification. This assures that 
all illuminated drops (contained within the light sheet) appear in the 
photographs with a sharp edge definition. The measurement volume is defined by 
the product of the projected frontal area of the nacelle and the light sheet 
thickness. Information on drop sizes and spatial distribution is also obtained 
from this photograph. 

When drops are photographed using the optical set up as shown in figure 1, with 
proper focusing and depth of field selection, the droplet edge definitions in 
the resulting photographic images are still far from ideal as shown in 
figure 2 . This is caused by two effects : 

1) The illumination of individual drops is not uniform, resulting in 
poor image definition of larger drops. 

2) Scattering and diffraction around smaller drops causes a "halo" 
effect around the drop boundary, causing the image to appear much 
larger than the actual size. 

An improved version of this technique utilizes laser induced fluorescence of a 
small quantity (less than 10 ppm) of dye (Rhodamine 6G or Fluorescene) 
introduced in the water. The fluorescence spectrum of the dye lies in a 
wavelength range longer than that of the incident light, therefore, the 
incident light scattered from the drops can be filtered out and only the 
fluoresced light is photographed. Shott colored glass filters are ideally 
suited for this purpose. This technique results in much better edge definition 
of droplets in the photographic images and a nearly uniform illumination of 
droplet interior (see figure 3) . 

An automated digital image processing system was developed to analyze the 
photographic images. The image processing algorithm detected drop edges, 

defined droplet boundaries, calculated the area of the drop images and an 
equivalent drop diameter based on this area. Droplet size distributions were 
constructed once statistics on a sufficient number of drops were available. 
The liquid water content was determined from the total volume of drops 
contained within the projected area of the nacelle, together with the knowledge 
of the incident light sheet thickness. A description of the image processing 
system is presented later in this report. 

For the determination of instantaneous mass flow rate of airborne liquid water, 
a measurement of drop velocities is needed in addition to the LWC measurement. 
This may be accomplished non- intrusively by laser double pulse photography of 
the moving droplet field. Two images of each drop appear in the photograph 
(see figure 4) and the drop velocity may be determined by the measurement of 
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Figure 1. Droplet Illumination and Photographic System 


- 4 - 



ORIGINAL PAGE IS 
OF POOR QUALITY 



Figure 2 . Droplet Photograph Without Fluorescent Dye in the Water 


- 5 - 



ORIGifML PASH 
OF POOR QUALITY 



Figure 3. Droplet Photograph With Fluorescent Dye in the Water 
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Figure 4 


Double Pulse Photograph of Droplets in a Moving Air Stream 



the distance translated in a known time interval. Solid-state lasers such as 
Ruby or Neodymium: yttrium aluminum garnet (Nd:YAG) lasers may be operated in 
the double pulse mode, wherein the. Q-switch is opened twice in rapid succession 
to split the energy of the flash tube between two intense pulses. The time 
interval between the two pulses may be adjusted in the range of 50 to 200 /is. 

The product of the LWC evaluated over the nacelle frontal area, the drop 
velocity and the nacelle frontal area is the instantaneous mass flow rate of 
airborne liquid water into the nacelle. 


WATER INGESTION SIMULATION FACILITY 

In the first phase of the program, an experimental facility was developed to 
simulate engine water ingestion and to calibrate the non- intrusive optical 
technique for the quantitative measurements of engine water ingestion. A 
description of the facility follows. 

The engine water ingestion simulation was set up in an open circuit wind tunnel 
of 18 -inch x 18 -inch test section. A schematic diagram of this facility is 
shown in figure 5. Air was supplied to the test section by a high capacity 
blower via a settling chamber, screens and a nozzle contraction. The blower 
output was adjusted by a damper vane ring at the inlet. A water spray nozzle 
was mounted on a streamlined sting support at the entrance to the test section. 
The nozzle consisted of seven 2 cm long tubes of 1.6 mm inside diameter. The 
water was supplied to the nozzle from a pressurized tank via a flowmeter and a 
ball valve. The individual tubes of the spray nozzle were adjusted such that 
most of the spray water entered the simulated nacelle with the tunnel running. 
The droplet-air mixture entering the simulated nacelle passed through a 
separator box which contained a series of baffle plates. The baffle plates 
subjected the flow to a series of sharp turns, thereby separating the droplets 
from the airstream and collecting them in a water film on the plate surfaces. 
Most of the water entering the nacelle in the form of droplets was collected at 
the bottom of the separator box. The air was evacuated from the separator box 
by means of the suction provided by the inlet of a secondary blower. The air 
from the exit side of the secondary blower was dumped back into the tunnel, 
downstream of the nacelle location. Calibration checks were made to determine 
the carry over loss of water in the form of fine drops. For the spray flow 
rate employed (300 cc/sec) , the carry over loss was about 5 percent, with 95 
percent of the water sprayed into the nacelle being collected at the bottom of 
the separator. With a steady spray, the total mass of water collected, divided 
by the time of collection was close to the time averaged mass flow rate of 
water entering the nacelle in the droplet -air mixture. 

The optical set up is shown in figure 1. The laser sheet of 9 mm thickness was 
produced in front of the nacelle inlet plane approximately 2 mm from the inlet 
plane to avoid illuminating the water film on the rounded nacelle entrance. 
The illuminated droplet field was viewed through the side of the tunnel. The 
camera axis formed an angle of approximately 20 degrees with the nacelle axis. 
In viewing the droplet field through the plexiglas wall of the tunnel at such a 
shallow angle, multiple reflections of illuminated drops were encountered 
within the plexiglas wall. To alleviate this problem, a window extension was 
mounted on the tunnel side wall at a 20 degree angle, such that the camera line 
of sight was normal to the window at the end of the extension. 
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Figure 5. Schematic Diagram of Engine Water Ingestion Simulation Facility 



The objective of the experiments was to compare the instantaneous water flow 
rate into the nacelle as determined by the present non- intrusive optical 
technique with the known time averaged water flow rate. Ten photographs were 
taken for a fixed spray rate and airspeed to yield a set of ten successive LWC 
measurements in a steady spray. The variation in the drop velocities within a 
given double exposed photograph or between two photographs at fixed tunnel 
speed and water injection pressure was found to be less than 5 percent. 
Therefore, a single drop velocity was used for water mass flow rate 
determination. A sample photograph of droplet field is shown in figure 6. The 
nacelle boundary was photographed separately and superimposed onto the droplet 
photograph in figure 6. Notice that the drops within the light sheet are 
uniformly illuminated and appear with sharp boundaries. Unfortunately, due to 
light scattered from drops contained in the light sheet and that from the beam 
dunfp and the transmitting side plexiglas window, some of the drops outside the 
laser sheet appear with faint images. During image processing of the photo- 
graphic negatives, these faint drops, which are present outside the light sheet 
may easily be eliminated by setting proper threshold criteria. 


IMAGE PROCESSING SYSTEM DESCRIPTION 

Processing system architecture is depicted in figure 7. Image acquisition, 
display and processing was accomplished using a De-Anza ID -5400 image 
processing system. The hardware package incorporates a vidicon and power 
supply for analog image formation, three image refresh random access memory 
channels, RAM, digital video array processor, and color video display. The 
analogue signal from the video camera can be digitized by an A/D converter and 
fed directly to the array processor which in turn controls the data flow and 
writes the data into one of the memory planes at a rate of 30 frames/sec. The 
digitization process converts each picture into 512*512 matrix element 
(pixels). Each pixel is one byte number (256 resolution level) representing 
the average optical density in an elementary cell, the size of which dictates 
the spatial resolution of the system. While digitization can proceed at video 
rates of 30 frames per second, a program is used which creates one digital 
frame from the average of 64 consecutive digitized frames . Thus , the image 
formed has a low level of random noise caused by the vidicon and the digitizer 
electronics. The digitized image information is then stored on a mass-storage 
device for further off-line processing. 

Software residing in the host computer (PDP 11/34) operates through a direct 
memory access (DMA), interface through which the PDP-11 sends and receives 
information from the video processor registers or from the RAM channels via a 
driver program. The vidicon image digitization, averaging and storage 
capability is part of this (DMA) interface software. 

To study a droplet picture , the negative is taped on to a transparency 
containing a 1 cm x 1 cm grid formed by fine lines. The grid lines are 
carefully oriented with respect to the droplet images on the negative such that 
none of the drop images is intersected by the grid lines. The negative, 
together with the grid lines, is then mounted on a flat light-table. The 
vidicon is focused on the plane of the background illuminated negative by means 
of a macro lens. The magnification on to the vidicon is adjusted such that a 
1 cm^ area of the negative enclosed by the superposed grid lines nearly fills a 
video frame containing the 512 x 512 pixel array. After accounting for the 
photographic magnification of the image on the negative and a specification of 
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Figure 6. Typical Photograph of the Droplet Field at the Nacelle inlet 
plane 
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Figure 7. Image Processing System Architecture 











4 pixels for the minimum size drop accepted, an overall resolution of 0.1 mm 
was estimated for the diameter of the smallest drop measured. 

The negative was scanned, one grid square at a time over the projected frontal 
area of the nacelle. The calculated drop volumes were then added up for all 
the squares to evaluate the liquid water content over the projected frontal 
area of the nacelle. This procedure ensured that every drop is accounted for 
and is counted only once. An example of a digital sub image representing one 
grid square is shown in figure 8. An enhanced image with defined droplet 
boundaries is shown in figure 9. 

Image processing software was developed to detect drop edges, define droplet 
boundaries , and calculate the area of drop images . The excellent contrast 
between the images of drops within the light sheet and the background allowed a 
simple thresholding criterion to define regions containing drop images. 

The various image processing programs along with computer program needed for 
droplet statistics determination are attached under Appendices A and B. 
respectively. ' 


RESULTS AND DISCUSSION 

Results are presented here for fixed values of spray flow rate and airspeed in 
the wheel spray simulation tunnel. The flow rate through the spray nozzle was 
maintained constant at 0.3 liters/sec, which resulted in a water jet velocity 
of 22 m/s based on the total flow area of the seven tubes in the spray nozzle. 
As discussed later, this flow rate was chosen to simulate extreme cases of 
engine water ingestion resulting from wheel spray. The airspeed was maintained 
at 61 m/s. The large difference between the speeds of the air and the water 
jets caused an aerodynamic breakup of the water jets into small drops, which 
accelerated along the flow direction to approach the airspeed. The measurement 
station was located 3 m downstream of the spray nozzle and immediately upstream 
of the simulated nacelle inlet. 

The droplet velocities were determined from a double exposed photograph of the 
moving droplet field as shown in figure 4. A 100 fis time interval between the 
two pulses was employed. The droplet velocities were found to be within 5 
percent for a large number of droplet pairs. Furthermore, within this small 
variation of droplet velocities , no correlation was found between the drop 
velocity and size. A single average value of velocity was assigned to all 
droplets. This average value was found to be only 41 m/s, i.e., 67 percent of 
the airspeed. Thus, the 3 m distance between the spray nozzle and the 
measurement station was insufficient to accelerate the drops to the tunnel 
airspeed. 

A series of ten photographs was taken during the period January -February 1986 
for the liquid water content determination. The width of the light sheet was 
maintained at 9 mm. The photographic negatives were analyzed by the image 
processing technique discussed in this report. The volume of a drop was 
calculated as the volume of a spherical drop having the same equivalent 
diameter. The equivalent diameter was calculated from the image area of each 
drop. For a non- spherical drop this procedure leads to a higher value of the 
calculated drop volume. The percentage error introduced therefore depends upon 
the how different the drop is as compared to the spherical shape. 
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Figure 8. Digital Subimage Representing One Grid Square 




Figure 9. Enhanced Version of Digital Subimage Shown in Figure 8 
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The results of the image processing analysis are shown in table I. The raw 
droplet data has been shown under Appendix C. For each of the ten droplet 
pictures analyzed, the total number of drops counted over the nacelle frontal 
area, the mean droplet diameter, the volume weighted mean diameter and the 
total liquid volume are presented. The standard deviation of the mean and 
volume weighted mean droplet diameters as calculated from the ten pictures was 
in the range of 7 to 8 percent of the average . This indicated that the 
picture -to -picture variation of the calculated mean diameter was ^latively 
small. The picture- to-picture variation of the total liquid volume was larger: 
the standard deviation of the ten measurements was 31 percent of the average. 
The average value of the total liquid volume from the ten pictures was used to 
calculate the average liquid water content, i.e., 


(LWC) a = = 450 gm/nr*; where ^w = 1.0 gm/cnr* 

Af . 6 

The average mass flow rate of water was then calculated as 

“w - (LWC) a . A f . V d 
“ 332.5 gm/sec 

The average mass flow rate of water in the moving droplet-air mixture as 
measured by the present non- intrusive optical technique thus indicates a flow 
rate approximately 10 percent higher than the actual value of 0.3 1/s. Reasons 
for this discrepancy between the measured and the actual mass flow rates are 
discussed below. It should be noted however, that there is a significant 
picture- to-picture variation in the instantaneous water mass flow rate 
determinations caused by measured variations in the instantaneous LWC values. 
The test droplet spray was quite dense in comparison with the certification 
requirement of 4 percent by weight of liquid water in the airstream. The 
present 0.3 1/s water ingestion rate in a 61 m/s airstream over the frontal 
area of the 15.24 cm diameter nacelle translates to a water flow rate/airflow 
rate ratio of 22.5 percent. Such a high water spray rate was used to simulate 
extreme cases of engine water ingestion resulting from wheel spray. 

There are three factors that affect the accuracy of LWC measurement by the 
present technique: 

1) Droplets largely outside the light sheet but at the boundary are 
partially grazed by the light sheet and show up on photographs, thus 
increasing the LWC measured. 

2) The procedure for calculation of drop volume from its non- spherical 
images on the photograph relies on an equivalent diameter which is 
calculated from the enclosed area of the image. This procedure tends 
to over-predict the volume of the non- spherical drop and hence leads 
to a higher value of measured LWC. 

3) The illuminated drops contained in the light sheet are viewed by the 
camera through a dense spray. Therefore, there is a possibility of 
some illuminated drops being masked by droplets present in the view 
path of the camera. This masking will result in a lower value of the 
measured LWC. 
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TABLE I 


Summary of Data Photographs 

V a — 61 m/s; V w - 21.7 m/s; q w = 300 cm 3 /sec; x = 3 m; Af - 182.4 cm^; 
5 = 9 mm 


Negative No. of Drops Mean Dia. Vol. Mean Dia. Total Liquid Vol. 


No. 

* 

mm 

mm 

mm ^ 

1 

267 

0.55 

0.89 

97.4 

2 

207 

0.58 

0.86 

68.2 

3 

256 

0.51 

0.87 

89.0 

4 

155 

0.58 

0.86 

52.3 

5 

179 

0.63 

0.92 

88.0 

6 

212 

0.56 

0.88 

75.0 

7 

243 

0.59 

0.91 

95.7 

8 

128 

0.50 

0.68 

21.5 

9 

337 

0.53 

0.80 

91.8 

10 

169 

0.60 

0.87 

57.8 



* Average: 0.56 

0.85 

73.7 


Standard Deviation: 0.04 

0.066 

23.0 

Std 

Dev. percent of Mean: 7.0 

7.8 

31.0 

NOTE: 

Mean Dia. = ^n^D^ . Vol. Mean Dia. 

SnjDj 4 



Sn^ 

SnjDi 3 



Total Liquid Vol. = Sn^ n D^ 3 




6 
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Despite these three sources of errors, the present non- intrusive technique is 
still considered a good method for measurement of the liquid water flow rate in 
a moving droplet-laden airstream. 

The droplet size distribution determined from the more than 2000 drops counted 
in the ten pictures is shown in figure 10. 


CONCLUDING REMARKS 

1) A non- intrusive optical technique has been developed for quantitative 
determination of instantaneous liquid water mass flow rate in a droplet 
laden airstream. The technique is generally applicable to the problem of 
water ingestion into an engine resulting in up to 22% water by weight 
(rain or wheel spray) . 

2) ' The technique yields instantaneous spatial distribution of droplets at the 

nacelle inlet plane as well as the droplet size distribution. 

3) Significant variation in the instantaneous values of the liquid water 
content was encountered in a droplet laden airstream produced by injection 
of a steady water spray in a steady airstream. The standard deviation 
among ten separate instantaneous determinations of the LWC was 31 percent 
of the average LWC. 

4) The average liquid water flow rate as determined from the average LWC and 
drop velocity measurements was approximately 10 percent higher than the 
actual spray flow rate. 
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FREQUENCY 17 .) 


I 



DIAMETER (MM) 

Figure 10. Droplet Size Distribution at the Nacelle Inlet Plane 
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APPENDIX A 


IMAGE ANALYSIS PROGRAMS 



grig3::al paslc 

OF POOR QUALITY 


FORTRAN IV-P V 02-51 10 : 58:17 26 -APR -35 P 1 

DROPS. FTN /TR : ULOCKS/WR 


0001 

C 

c 

c 

0002 

0003 

0004 
*0005 
0006 
0007 
00011 

0009 

0010 
0.0 1 1 
0.0 1 2 
0013 
00 1 4 
0015 
00 1 6 

00 1 7 
0 . 01 H 
0.010 
0020 

002 1 
0022 

0023 

0024 

0025 

0026 

C 

C 

C 

C 

0027 

0025 
0023 
0030 
003 1 

0032 

0033 

0034 

0035 

0036 

0037 
0036 
003 !) 

0040 

004 1 900 

0042 

0043 


PROGRAM DROPS 

Drops analysis based on averaged maximum gradient 

LOGICAL*! ANSA ,AN30 
•JYTI; 2IIU0 

REAL* i) II I ST! 0:255 > ,DMAXO 
BYTE HUP IC! 5 12 ) , RtJGRA< 512) 

BYTE PIXEL! 5 12, 16 ) 

COMMON /P IMTA/P IX EL 
IHTEGE R*2 IVLT<256> 

I NTHGE R*2 ARM IN 

I HTE GER*2 IBUEOf B ) 1 Basic data processing parameters 
EQUIVALENCE ! IHUEO! 1 ) , KREDC ) 

EQUIVALENCE < I HUFO< 2 ) , SCALE ) 

EQUIVALENCE ! IHUEO! 4 ) .ARM 1 N ) 

EQUIVALENCE ! I BUEO! 5 ) , ICORTA ) 

INTEGCU*2 SL , S3 , HL , NS 
REAL. *4 A HUE < 200 ) , P HUE ( 200 ) 

INTEGER *2 LIMIT! 4,200) 

COMMON /M 1 NARE /AM I NAR 
INTEGER* 2 IHUEE! 3 ) , PL , PC , CL , CC 
REAL *4 AREA, PER IM 
EQUIVALENCE < I HU EE ! I ) ,PL ) 

EQUIVALENCE ! IHUEE!2),PC) 

EQUIVALENCE < IHUEE! 3 > ,CL ) 

EQUIVALENCE < IHUEE! 4 >,CC) 

EQUIVALENCE ! IHUEE! 5 ) .AREA ) 

EQUIVALENCE < IHUEE < 7 > , PER IM ) 


Open data file 

OPEN ! UN IT= 3 , NAME 3 * SY : C 300 , 300 ] RAIN. DAT ' .TYPE- ‘OLD’ , 

1 I) ISP 3 1 KEEP ' .ACCESS* * DIRECT 1 , RECORDS I ZE = 4 , FORM* * UNFORMATTED 1 ) 

READ! 3' 1 ) IBUFU 

IF ! KREDC, HE , 0 ) GO TO 900 

KREI)C= 1 J Initial record for parameters storage.no data 
ARM IN -51 Default min area 
I CURTA -6 
SCALE * 1 

TYPE * , ' Processing parameters Initialization, defaults are:' 
TYPE *,' Min area 5 Sq. pixels. Scale 1 mm/pixel ; Change <Y/N)I* 
ACCEPT 101 .ANSA 
IE (ANSA. HE. * Y * ) GO TO 111 

TYPE *,’ Enter Min area !Sq, plx) and Scale {min/plx)' 

ACCEPT *,ARMIN, SCALE 
GO TO 111 

TYPE File has already data. Stops execution !Y/N)' 

ACCEPT 101 , ANSA 

IE < ANSA , EU . *Y ’ ) GO TO 950 


C 

0 First detector 

C 

0044 111 TYPE *,' Begin coarse detection* 
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ORIGINAL PAGE 
OF POOR QUALITY 


FORTRAN IV-F J V/tt-M li7: GO : 17 2G-APR-05 i 2 

DROPS. FTN / TR s ULOCKS/VR 


niur> 

OlMi) 

•ihmu 

dim*) 

0000 

ooo \ 4 

0002 j 

0003 

c 

c 

c 

0004 s 
0000 


noor, 

/ 

00 OU 

oooo 

0000 

O0G1 


O0G2 

ouc>3 

ooGi 

0000 


noor, 

000 7 


0"GU 
0000 
00 /ir 

007 t 

on/;? 

007 1 
0074 


007 r > 
0070 
00/7 
0070 
0070 
OOSlO 
non ( 
oo!j;» 
00ii:i 
00!) 4 


10 

c 

c 

c 


c 

c 

c 

c 

c 

c 


20 

i; 

t: 

c 


r,o 


*31 


CAM. SETUP 

TYPE * , ' Smoothing before compute gradient (Y/N)' 
ACCEPT 101 .ANSA 
II* < ANSA. El). , H ) GO TO 4 
CALL SMOOTH 
ISO TO *3 
1)0 :) IC»l ,20G 
1VLT( K )«K/2 
CALL VL TCOIU 1 . IVLT ) 

Save original picture 


CALL SETUP 

OPEN < Ull IT* 7, TYPE® * SCRATCH 1 .ACCESS® 'DIRECT* , I N IT I ALS IZE*5 1 2 , 
1REC0RUS12E-12U) 

MI)LK»1 

DO 10 K-1,32 
1. IN *<K-1 >*10+1 

I COD® IMAGE < 0, 1. IN . 1G . 'R ' .P IXEL ) 

ISTAF®.)ESC{ 7. P IXEL , NBLK . 409G , 1 < NBLK ) 

CUNT I HUE 

Transfer smooth picture to channel 0 


CALL SETUP 
CALL I MSEL (2*0,0) 

CALL OUTHIILU) 

CALL KHOPER 

Compute gradient 

CALL SETUP 
CALL GRAD 1 E 

Restore original picture 

CALL SETUP 
HULK"! 

DO 20 K :s 1 ^ 32 
L !H«f K-l > MG+l 

ISTAT«*JLEEf 7 , P IXEL .NBLK . 4096, 1 .NBLK) 

1 COD* l MAGE < 0 . L 1 H . 1 G • V .PIXEL ) 

CONTINUE 

Threshold gradient picture 

CALL SETUP 
CALL IHSELU.B.0) 

CALL OUTHHI. { 4 > 
i:ALL KHOPER 
CALL SETUP 

TYPE Threshold gradient Image '. ICORTA 

CALL SETUP 

DO 01 JI.T» 1 « ICORTA 

IVLT 1 Jl.T )«0 

DO 02 JLT» ICORTA* 1 . 20G 
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ORIGINAL PAGE & 
OF POOR QUALITY 


FORTRAN 

IV- 

S V01-51 10! 50 ? 17 2G-APR-05 

DROPS. FTM 

/TR 2 ULOCKS/WR 

0005 

SI 

I VI. 7 < »ii.t >»isr, 

0006 


CAM. V IDTIIRI 1 , 1VI.7 ) 

i)iMi / 


TYPE Accept throshold <y/n)’ 

DDini 


ACCEPT 101 ,ANSA 

iron*) 


IK (ANSA. HE. Ml’ ) CO TO G0 

0000 


TYPE Enter gradient Imago threshold' 

!)i)') t 


ACCEPT * , IC0R7A 

))')•):! 


CO TO 50 

Di)*)! 

101 

KOKH AT ( A 1 ) 

1)1)* >4 

S0 

CALL SETUP 

oins 

* 

DO 70 JI.T-* IC0R7A+ 1 ,256 

OD'lij 

70 

1 VLT < *)LT MJl.T-l 

on*)/ 


CALL VLTCOIU 1,1VLT) 

00’JO 

C 

i; 

c 

CALL SETUP 


Save gradient picture 

01)*)*) 

NISLKM 

1)11)0 


DO 1)0 K= 1 , 32 

own 


1. 1M«< K- 1 >*16+1 

did;! 


I COD 5 * IHACEI 1 , L I M . 16, *R\ PIXEL) 

010! 


l STAT*3ESC{ 7, PI XEL , MBLK , 4056 , 1 , f|BLK ) 

0 104 

00 

CONTINUE 


c 

c 

Remove background variations 

0105 

TYPE Remove background variations <Y/N>' 

0 I0ii 


ACCEPT 10l,AH5U 

010/ 


IK {AHSO.ME.MD CO TO 05 

0100 


CALL I USE l.< 2,0 ,0) 

010’) 


CALL OUT Hill. < 4 ) 

01 10 


CALL KHOPER 

Dili 


CALL SETUP 

0 1 1 1 


CALL IHSEL< 1,0,0) 

011! 


CALL OUTNI1LI 2 ) 

Dili 


CALL KHOPER 

0 1 IS > 


CALL SETUP 

0 l Hi 


co ro do 

0 l 1 7 

os 

I SI)™ 100 1 Variations scale 

Dill) 


CALL IIICPASI ISO, 1 ,0) 

urn 

r: 

c 

c 

CALL SETUP 


Restore gradient picture 

0 l .10 

HI) LIC n 1 

0 1 1 l 


DO 120 K-1,32 

0111 


L III»<K-1 

0111 


l STA T«dLEE< 7 f P I XEL , HOLK , 4056 , 1 , NBLK ) 

0111 


I COD** IMAGE ( 2 , L I H , 1 G , ’ V ' , P I XEL ) 

0 1 IS 

110 

i % 

CONTINUE 


i; 

i* 

Compute Histogram 

011 S 

i if) 

2ERO-*0 

0 117 


TYPE * , ' Computing modified PDF' 

0111 ) 


DO 00 K =0,255 
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10150:17 


26-APR-OS 


FORTRAN IV- S V/T^-M 
DROPS. FTN / riCULOCKS/WR 


£ 4 


0120 

ou 

11 1 ST < K )»0. 

0 i 30 


Ill) 150 K™ 10,500 

1) 131 


n:ni >» ihaishc mm.'R 1 ,nupic> 

0132 


n:uti»iHAi.i:(2,K. i, nr ,uugra> 

m 33 


1)0 1 50 .1- 10, 5011 

m.n 


IF <liUGRA<OI.EU.2ERO) GO TO 150 

in 35 


A~IV< lUJGRAI J 1 ) 

0135 


! l"IV<UUPIC< J ) ) 

il 1 3 7 


11 1 ST < 1 I MllSfl 1 1 KI I*A 

1) iju 

150 

i: 

i; 

C 

CONTINUE 


D lap lay histogram and threshold determination 

nm 

l)HAXO*0. 

m in 


DO HID K«0,255 

III 4 1 


IF < II l ST { K I .GT . DMAXO ) DMAX0*H 1ST < K ) 

dm;! 

100 

CONTINUE 

Il l i l 


1)0 1 U2 K rt 0, 255 

n i 4 4 


II 1 ST < K ) "H 1 ST ( K ) /DMAXO 

II 1 4?» 

102 

CONTINUE 

0 l 4 5 


CALL IHSELi 1 ,U,0> 

11147 


CALL OUTNIIM 4 ) 

iimh 


CALL FHUPEK 

III 4*1 


CALL SETUP 

I) l fill 


1)0 150 IC« 1 • 25C 

i) i 5 i 

1G0 

IVLT(K)=0 

o i 


CALL VIDTIIIU 1 , IVLT ) 

.0 l •> :l 


If 1111=0 

1) 154 


CALL IMSIIISMIIST, 1THR) * 

i) l *i!i 


TYPE * , ’ Enter (1) 3lmp1e, {2) bl modal distribution 

0 155 


ACCEPT *. INODES 

ill 5 / 


CALL SEUISOtlflST. IT1IU , IMODES > 

1) 150 

300 

TYPE * • ’ Estimated thresho Id * « ITIIR 

1) l 50 


CALL SETUP 

II l till 


DO 301 JI.T-1, ITIIR 

II l 5 l 

301 

IVLT < %)LT )»0 

Dili;! 


DO 302 »)LT a ITHU+ 1 , 256 

1) 153 

302 

IVLTl JLT >«255 

II Hi 4 


CALL VLTCORM.IVLT) 

0155 


CALL SETUP 

0155 


IVI.TI l I^ll 

Dili 7 


I VLT ( 255 )** 1 

01 5 U 


CALL VLTCOIK 1 , IVLT ) 

01 50 


CALL SETUP 

D 1 7 il 


CALL l!URDEA< IVLT ) 

D l 7 l 


CALL SETUP 

nr/:! 


CALL D 1SII IS< II I ST * ITHR ) 

il i 7 1 


TYPE * , ’ Accept threshold (Y/N)' 

11174 


ACCEPT 101 ♦ ANSA 

0175 


IF (ANSA .HE . 'll’ ) GO TO 350 

D l 7 1» 


TYPE * « ’ Enter threshold 1 

l> l 77 


ACCEPT Mfllll 

D l 711 


IE { AIISO. HE . *ir 1 GO TO 320 

Il l 7*1 


CALL IHSELI 1 ,H«0> 

III HD 


CALL OOTHIIL < 2 ) 

nun 


CALL FHUPLR 


ORIGINAL PAGE fS 
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FORTRAN 

IV- 

S V02-O1 10:00:17 2G-ArR-05 P 


5 

DROPS. FTM 

/ TR:UL0CK3/VR 



II UK'! 


CALI. SETUP 



JIUII 


<•0 to jnn 



hum 

320 

CALL IIIi;PASUSOMt0) 



OHIO 


CALL SETUP 



nuui 


GO TO 300 



iiui; 

350 

TYPE *,* Default region the entire picture* change 

< Y/N ) 

in nil 


ACCEPT 101.AHSA 



nun 


IK < ANSA . EO . ’ Y * ) GO TO 400 



ll l 'HI 


si.**i 



nun 


30 s * 1 



IIUK! 


HI. =*010 



null 


ns=oiu 



HUM 


GO TO 400 



n i'is 

400 

TYPE *,* One curaor to define rcglon/pusli white key* 


n i 'jo 

i» 

CALL KEG I0< SL , 00 . HL , NS ) 



on 7 

i, 

400 

AM IHAR^ARH IN 



nui« 


TYPE *.* Region definition SL-'.SL.* SS-^SS.' HL = 

4 

ML. 1 

nui'i 


I VLT < 1 )»0 



0200 


I VLT < 2 )=*200 



n; !o 1 


I VLT { 200 ) lt 205 



n;!.n:! 


CALL VLTCOIU 1 , I VLT ) 



n 2 . 0:1 


CALL SETUP 



n;!iM 


TYPE * • * Computing drops size' 



0200 


CALL GOTA3< 1 , SL , 00 , HL , NS ,ADUF , POUF , L IM IT « NUMOOO ) 



n:!nij 


I VLT { 1 )~0 



0207 


IVLT< 20b >«1 



02.00 


CALL VI. TCOIU 1 , 1VI.T ) 



n;:.0'i 


CALL IIORUEAUVLT) 



IK! 10 


CALL SETUP 



IK! 1 1 


TYPE * , ' End coarse detector* 



0212 

000 

i' 

CLOSE { UN IT *7 ) 




c 

p 

Save ro3ult3 on dl3k 



02 l :i 

2000 

HUMOR *0 



|K!U 


3CALE2»3CALE*3CALE 



IK! 10 


TYPE *,* DATABASE: Initial o lenient \ KRE DC 



IK! 1 it 


DO 2001 K=* 1 , NUMUO J 



IK! 17 


-ItUHIIIMIUMUR+l 



IK! HI 


KUEDOICUE l)C> l 



IK! Ul 


PL«LIHIT< 1 , K ) 



0220 


PC-L I M I T < 2 , K ) 



ik!;! 1 


CL w L I H IT { 3 t K 1 



ll ;?;*:! 


CC«L IM ITC 4 .K » 



n:!;!:i 


AREA^ARUKI K SCALE 2 



n;!;! I 


PER IM^PIIUF ( K I* SCALE 



1)220 


WRITE Cl KRE DC ) IBUKE 



n;!:!0 


IF (ICREUC.EQ. 1201 ) GO TO 2025 



IK! 2 7 

2001 

CONTINUE 



ikmh 

2020 

TYPE *.* DATABASE : Final element * ,KREDC-1 



0220 


TYPE Number of objects In this f rawo ' , NUMDR 



ik! in 


W’R!TI;<3*1) I I1UFO J Update general Information record 



l 

300 

CLOSE! UNIT® 3 1 



ik:k! 


EHD 




ORIGINAL PAGE ES 
OF POOR QUALITY 



FORTRAN * IV-PL 
GOTAS.FTN 1 


V02-51 

/TR : BLOCKS/VR 


14109:02 


25-APR-85 


SUBROUTINE GOTAS< CHAN, SL , SS , NL • NS , ABUF , PBUF , L IH IT . NUMOBJ > 


Drops geometrical caractor 1 zat 1 on 

INTEGERS CHAM, SL.SS.NL, NS 
BYTE I' IXEL<5!2, 15) 

COMMON /PINTA/PIXEL 
I NT EGER *2 FLAG( 5 12 , 2 ) . ID<30> 

INTEGER *2 ISL3( 90) , IS3S< 90 ) , EL ( 90 ) , ES ( 90 ) 
INTEGER*? CMFIK 90.2 ) , BEG IN< 90, 2 ) , I END < 90, 2 > 
I NTEGER*2 OLD 1 1), NEW ID 


0009 


INTEGER* 4 HP 1X3(90) 

00 10 


K EAL *4 K P E R < 90 ) 

00 I 1 


REAL *4 AI1UF < 1 ) , IM1UT < 1 ) 

0012 

p 

INTEGER*? L I H I T < 4, 1 ) 


C 

p 

Initial 1 za 

0011 

V* 

CALL 2 !A< FLAG( t, 1), 1024) 

00 14 


CALL 2IA(RPER, IDO) 

0019 


CALL 7. 1A< HP IXS , 100 ) 

00 Ui 


CALL 2 1A< EL , 9 If ) 

»w t / 


CALL 2 I A ( E 1 , 90 ) 

0010 


CALL 2 I A( CNTIU l, 1 >, 100 > 

00 19 


CALL ZIA<UE<5IN< 1,1), 1U0) 

0020 


CALL 2 IA( IE HIM 1 , 1 ) , 1 00 ) 

002 l 


CALL I T I A < 12000, I SI. 3, 90) 

0O22 


CALL IT I A ( 3 2000 , I S SS , 90 ) 

002 1 


MOHORJ-W 

0024 


I :»W“ 1 

0029 


NEWC-0 

002G 


1101. C *0 

0027 


01.1)1 1)"! 

0020 


H3S-S3+N3 

0029 


II3I.-SI.-1 

0030 

p 

HLL-H3L+NL 


i; 

c 

Main loop 

001 1 

1)0 1WW L- 1 ,HL, 16 

0012 


lltl.OCIC-M IM0< 1G.ML-L ) 

0013 


1. IHE-l+MSI. 

0014 


|i:0l)«IH/\0E<CII/\H,LIME, IBLOCK, 'R' .PIXEL > 

0019 


1)0 100 KK“ 1,10 

0010 


H1-3S 

0017 


1. IIIM-L1NE ►KK-1 

0010 


MI-VOW 1 Start searching for an object 

0019 

ms 

110 11W .1-N1.M33 1 tool: for left edge 

00 40 


II- < PIXEL! J,KK>) 00 TO 120 

004 1 

l UJ 

CONTI HOC 


c 

Mo wore edges on this lino 

004 2 


00 TO Uli 

0041 

120 

lll=.J 1 Left edge 

004 4 


DO 1J1 .1 = 111, H33 1 Look rlgh edge 

0049 


IE ! .HOT. PIXEL! J.KK ) ) 00 TO 14W 

00 4 IS 

130 

CONTINUE 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


FORTRAN IV-Pl V02-51 14:09:02 25-APR-05 PA 12 

GOTAS •PTN TK : ULOCK3/VR 


0047 

0040 140 

i; 

0 

c 

0049 

0050 

0051 

C 

C 

0052 

0053 

0054 101 
C 

C 

0055 

C 

C 

0055 

i; 

i: 

0057 102 

00511 

005*1 

0050 

005 1 

C 

C 

0052 

0053 

0054 

0055 
0055 
0057 

C 

C 

0050 

0050 


C 

0070 103 


007 1 
0072 
007 3 

007 1 100 

00/5 

00/5 


•l“M33+l J High edge Is the picture boundary 
H2-J-1 


Kind ID of this object from previous flag lino 

CALL MATCH! FLAG! 1 , 3- 1 3V ) , N 1 , N2 . HUM , I U ♦ HS3 ) 
m : .wii)«iD( 1 ) 

IK < HUH • LT . 2 ) GO TO 101 

If many II) condense them Into one 
CALL COM ! D( UPEIt , HP I X3 , 1 3L3 , 1333 , EL ,E3 • HUM , ID ) 

«;«) TO 102 

IF ( HUM . ML .0 ) GO TO 102 

If thro re Is not ID them creates one 
CALL K IHD( 90, HEW ID , HP I X3 , ULU 1 1) ) 

begin count Ing perimeter new object 
RI’Elt! NLWID ) a H2-N 1+2 

And update flag buffer and counters 
CALL IT IA( HEWI D • FLAG! M 1 , I3V),H2-Ml+d ) 

ni:wi>i<i:vc+i 

CM TR< HI- WC , I3V>«M!:VID 
DEC IM< HEWC . ISWKIU 
I EMU! HE VC ♦ I SW >«H2 

Update statistics buffer 

HP IX3< MEWIO )*»HP I XS! NEWI U )+N2-N 1 4 1 

IF i l5LS!MKVIi> >.Gf J.INIA > 131. 5 ( NEVI D >*L I H I A 

IF < I333M4EWIU) .GT.H1 ) I S5S< NEWI D >**H 1 

IF M:3<M:VII)).LI\N2> ES!NLVID)*H2 

i:i.( HLWII) )-L I H I A 

CALL PER IN! NOLC # MFWC «CHTR! 1 « I5W > , HEV1 D , RPER , DEG IN! 1 ,3-ISV>« 

1 BEGIN! 1,1 3W ) , IEND< 1 , 3- I SW ) , 1EMD! 1 , I3W) ) 

Send bach for more data on the same I I no 
H 1 =M2 + 2 

IF ! M 1 . LE * NS5 ) GO TO 105 
Search buffers for terminated drops 

CALL EIIOFL I ! RPER , HP 1X3 ♦ ISL3 , 1333 ,EL ,E3 , NLL , L IH IA,CMTR< l , 3- ISV> , 

1 Cl« f R< 1, 1 SV ) , NOLC , NEWC , NUMOU J , DEG I N f 1 , 3- 1 3V > , DEG I Ml 1 , I5V>« 

2 I END! 1 ,3-ISW), I END! I , 1 3W ) ,AWJF , PIJIJF , L IM I T ) 

CALL Z I A! FLAG! 1 , 3- ISW) ,5 12 > 1 Zero flag line 
ISW b 3-ISV 1 Switch the flag 

IH)M>Mi:WC 

CUHT (HUE 1 Cto3e main loop 

RETURN 

EMU 
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19101:30 


04-APR-05 


FORTRAN ’IV-F V02-51 

HATCH . FTN /TR s BLOCKS/VR 


0001 


0002 

0003 

0004 

0005 

0006 
0007 
0000 
0009 


SUBROUTINE MATCH < FLAG , N 1 , N2.NUM, ID.NS ) 

C 

C RETURN 1*0. OF OBJECT FOUND ON PREVIOUS LINE 
C FLAG - OLD FLAG BUFFER 
C N 1 - BEGIN, N2 • END OF OBJECT 
C NUN - NUMBER OF OBJECTS LOCATED 
C ID - I . D . OF OBJECTS LOCATED 
INTEGER FLAG(I) 

INTEGERS ID< 1 ) 

M 1 * N i - 1 
M2 - N2 ♦ 1 
IF (Ml.LT.l) Ml - i 
IF < M2 . GT . NS ) M2-NS 
NUM - 0 
K * 0 


C 

C FIND ALL FLAGS 
C 


0010 


DO 

10 J-Ml, 

M2 





0011 


IF 

<FLAG< J ) 

.EQ 

.K> 

GO 

TO 

1* 

0012 


IF 

< FLAG! J ) 

.EQ 

. 0 ) 

GO 

TO 

10 

0013 


K - 

FLAG( J ) 






0014 


NUM 

- MUM ♦ 

1 





0015 


IIK 

NUM) - K 






0016 

10 

CONTINUE 






0017 


IF 

(HUM.LT. 

2) 

RETURN 




C 

C REJECT DUPLICATES 


C 


0018 


CALL CORSRT< ID, ID, NUM ) 

0019 


N - 1 

0020 


K - ID( 1 ) 

002 1 


DO 20 J«2 , NUM 

0022 


IF < K . EQ . ID< J ) ) GO TO 20 

0023 


N - M + 1 

0024 


IIXM) • ID< J ) 

0025 

20 

CONTINUE 

0026 


NUM - N 

0027 


RETURN 



ORIGINAL PAGE IS 
OF POOR QUALITY 


FORTRAN IV-r 
PER JM.FTN 


V/J2-51 19101 >04 04-APR-8B 

/TR:BLOCKS/VR 


P 1 


0001 


0002 

0003 

0004 


SUBROUTINE PER IM< NOLC .NEWC.NECNTR, NEWID.RPER .OBEG IN .NBEG IN. 

1 OEND.NEND) 

C 

C THIS SUBROUTINE ACCUMULATES PERIMETER MEASUREMENTS FOR THE ’OLD’ 
C LINE BETWEEN N1 AND N2 OF THE 'NEW* LINE. 

C 

INTEGER NECNTR! 1 l.OBEGIN! 1 ) , NBEG IN( 1 ),OEND! 1 ) ,NEND( 1 ) 

REALM RPER( 1 > 

IF < NOLC . EQ .0 ) RETURN 
C 


C SEARCH 'OLD ' LINE FOR MATCHING OBJECT SEGMENTS. 
C 


0005 

0006 
0007 

0000 

0009 

0010 
0011 

0012 


LI - 0 

DO 100 J-l.NOLC 

IF (UEND! J l.GE. NBEG IN{ NEWC 1-1. AND. 

1 OBEGINIJ > . LE . NEND< NEWC ) -*■ 1 > GO TO 101 
GO TO 100 

101 ONEWC - NEWC - 1 

IF < ONEWC.LT. 1 ) GO TO 301 

IF < NEWID . EQ . HECNTRI ONEWC > .AND .0ENDIL2 l.GE. NBEG I N< NEWC ) 1 
1 GO TO 302 
GO TO 301 
C 

C IF PERIMETER HAS BEEN ADDED FOR PREVIOUS PARTICLE SEGMENT ON 
C HEW' LINE WITH SAME I.D. AND THE 'OLD' MATCHING SEGMENT 

C OVERLAPS BOTH 'NEW SEGMENTS. THEN THE OVERLAP IN THE 

C CURRENT SEGMENT IS SUBTRACTED FROM THE PERIMETER VALUE. 

C 


0013 

302 

** lT> 

il 

301 

0016 

303 

0017 

0016 

0019 

100 

0020 

0021 

0022 

0023 

0024 

102 

0025 

0026 

104 

0027 

0026 

103 

0029 

0030 

200 

0031 

0032 



RPER! NEWID ) - RPER! NEWID 1 - SORT { < OEND( L2 1 

1 - FLOAT! HE HD (ONEWC >11**2 + 1.) + NBEG IN! NEWC) - HEND! ONEWC) 

GO TO 303 

RPER! NEWID ) ■ RPER! NEWID ) + SORT! ! NBEG IN! NEWC ) 

1 - FLOAT ! OBEG IN! J > ) )*"2 ♦ 1 . ) 

LI - 0 
K • J ♦ 1 
GO TO 102 
CONTINUE 

IF ! L 1 .EQ.0) RETURN 
L2 - 0 

IF ( K .CT .NOLC > GO TO 200 
DO 103 J-K.NOLC 

IF !OEND! J ).GE. NBEG IN! NEWC )-l. AND. 

I OBEG IN! J ) .LE .NEND! NEWC )+ 1 ) GO TO 104 

GO TO 103 

RPER! NEWID > • RPER! NEWID ) + (OBEGIN(J) - OENDIJ-l)) 

L2 - J 
CONTINUE 

IF ! L2 . EQ .0) L2 - LI 

RPERIHEWID) - RPER! NEWID ) ♦ SORT! ! NEND! NEWC ) 

1 - FL0AT!0END(L2)>>-*2 ♦ 1.) 

RETURN 

END 


-A9- 



FORTRAN' IV- 5 VJF2-51 18 » 59 t 52 04-AFR-85 t 1 

CONID.FTN /TR « BLOCKS/VR 

0001 SUBROUTINE CONIOI RPER ,NP IXS, ISLS, ISSS.EL .ES.NUH, ID ) 

C 

C USED TO CONCATINATE MANY I.D.’S 
C 


00*72 


INTEGER OLDID , ISLS< 1MSSS1 1>,EL< 1 ).ES( 1 ) 


0003 


INTEGERS NR I XS( 1 >♦ ID< 1 1 


0004 


REALM RPER< 1 ) 


0005 


NEVID - 10(1) 


0006 


DO 10 J*2 , NUN 


0007 


OLDID - ID<J) 


0008 


RPER< NEVID ) - RPER< NEVID ) ♦ RPER< OLD ID )• 


0009 


RPER< OLDID) - 0. 


0010 


NP 1XS( NEW ID ) • NPIXS(NEWID) 4 NPIXS(OLDID) 


0011 


NP 1 XS< OLD 1 D ) - 0 


0012 


IF < ISLS<NEVIU).GT. ISL3<0LDID>> ISLS(NEWID) - 

ISLSIOLDID) 

0013 


ISLS< OLDID ) • 32000 


0014 


IF < ISSSINEVIDKGT. ISS3<0LDID> > ISSSINEVID) - 

I SSS( OLD ID ) 

0015 


ISSS(OLDID) • 32000 


0016 


IF ( EL < NEW ID ).LT. ELI OLDID)) EL< NEVID ) - ELIOLDID) 

0017 


ELIOLDID > - 0 


0010 


IF < ESI NEW ID ).LT.ES< OLDID > ) ES(NEVID) - ES(OLDID) 

0019 

10 

ES< OLD ID ) ■ 0 


0020 


RETURN 


0021 


END 
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PORTRAN*IV-P 

V02-51 19:02:52 04-APR-85 

P 

i 

PCNCAT 

.FTN 

/TR : BLOCKS/VR 


000 1 

r 

SUBROUTINE PCNCATI NCEN.PCNTR, BEGIN, END. NEVID, 

RPER ) 



C THIS SUBROUTINE CONCATINATES ALL ENDING PERIMETER 

VALUES 

FOR EACH 


c 

r 

I.D. ON AN ENDING LINE. 


0002 

v« 

INTEGER PCNTR( 1 ), BEGIN* I ) , END< 1 1 



0003 


REALM RPER< 1 > 



0004 


DO \00 0*1 , NCEN 



0005 


IE < NEVID . EQ. PCNTR( 0 ) ) GO TO 101 



0006 


GO TO 100 



0007 

INI 

NPER< NEWID 1 - RPER< NEVID ) ♦ END< 0 1 - BEGIN(O) 

♦ 2. 


0000 

10/J 

CONTINUE 



0009 


RETURN 



0010 


END 




-All- 



19:03i*O 


04-APR-B5 


P 


FORTRAN 1V-P 

FIND.FTN 


V02-51 

/TR : BLOCKS/WR 


1 



0001 


0002 

0003 

0004 

0005 
000ft 
0007 
0000 

0009 

0010 
00 1 1 
00 1 1 ! 
0013 
00 1 4 
0015 


SUBROUTINE F IND( NBIN«NEWID,NP IXS .OLD ID ) 
C 

C USEO TO FIND A NEW BIN POSITION 
C NOIN - LENGTH OF BIN BUFFERS 
C MEWID - NEW BIN VALUE RETURNED 

C NPIX3 - SUN OF PIXELS BUFFER 

C OLDID - LAST FOUND I.D. 

INTEGER OLDID 
INTEGERM NPIXS< 1 ) 

DO 100 0-OLD ID , NO IN 
IF <NPIXS<0 1.EQ.0) GO TO 200 
100 CONTINUE 

DO 150 0-1, OLDID 
IF ( NP I XS < J > . CQ.0 ) GO TO 200 
150 CONTINUE 

TYPE *, • All bins filled' 

STOP 

200 NEW ID-0 

OLD ID-0 
RETURN 
END 
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FORTRAN '|V*-P V02-51 19:03:28 04-APR-85 P 1 

CORSRJ.FTN /TR :BLOCKS/VR 

0001 SUBROUTINE CORSRT< KEY.PTR.LEN > 

C 

C IN CORE SORT ROUTINE. SORTS THE VECTORS KEY AND PTR INTO ASCENDING ORDER OF 
C KEY. THE VALUES IN KEY ARE TREATED AS LOGICAL QUANTITIES HENCE SIGN BITS 
C HU3T BE TREATED ACCORDINGLY. ALPHABETIC RECORDS HAY BE SORTED SINCE 
C CHARACTER CODES ARE IN LOGICAL ORDER. THE PTR ARRAY CAN THEN BE USED TO MOV 
C RECORDS IN A DISK FILE. 


0002 

C 

IMPLICIT INTEGER! A-Z > 

00153 


IHTEGERM IMSK , JMSK , SMI32>,PTR< 1 >,TEMP 

0004 


LOGICAL *4 IMSL ,JMSL ,KEY! 1 1 , LTEMP 

0005 


DIMENSION HDRY(3,32),CBDI32) 

0005 


EQUIVALENCE I IMSK, IMSL 1 , ( JMSK , JMSL ) 

0007 


DATA SM< 1 1 , SMI 2 1 , SMI 17 1/020000000000,010000000000,0100000/ 

0000 


DO 1 1-3,16 

0000 

1 

SMI U - SMI 1-1 1/2 

00 10 


DO 2 1-10,32 

00 l 1 

2 

SMI I ) - SMI i-1 >/2 

0012 


LEV - 1 

0013 


BDRYI 2 , LEV 1 - 1 

0014 


BDRYI 3 , LEV 1 - LEM 

00 1 5 


CUD! LEV 1 - 2 

00 IG 

72 

CU - CHOI LEV ) 

00 1 7 


PL - BDRYI CH ,LEV ) 

00 1 0 


PU - BDKYICB+ 1 , LEV 1 

0019 


IK (PL.GE.P0) GO TO 75 

0020 


IMSK • SMI LEV ) 

0021 

81 

JMSL • KEY! PL 1 .AND . IMSL 

0022 


IK I JMSK , HE ,0 ) GO TO 05 

0023 


PL - PL ♦ I 

0024 


IK IPL.EO.PU) GO TO 73 

0025 


GO TO 81 

00 2 G 

85 

JMSL - KEY! P0>. AND, IMSL 

002/ 


IK I JMSK . HE ,0 1 GO TO 86 

002O 


LTEMP - KEYIPL 1 

0029 


TEMP - PTRIPL) 

0030 


KEYIPL 1 - KEYIPUl 

0031 


PTRIPL ) - PTIUPUl 

0032 


KEYIPU) - LTEMP 

0033 


PTRIPU) « TEMP 

0034 


GO TO 81 

0035 

8G 

PIJ - PU - 1 

003G 


IK IFU.NE.PL) GO TO 05 

0037 

73 

JMSL - KEYIPL ). AND. IMSL 

0030 


IF (JMSK. HE. 0) PL - PL - I 

0039 


IK < LEV .GE . 32 1 GO TO 75 

0040 


CH - CUD! LEV) 

004 1 


LEV • LEV ♦ I 

0042 


i;hii< lev i - i 

0043 


IJDKY I 1 , LEV ) - BDRYI CD , LEV-1 1 

0044 


HDRY I 2 , LEV > - PL 

0045 


BDRYI 3 , LEV 1 - BDRYI CB* 1 , LEV- 1) 

004 G 


GO TO 72 

0047 

75 

IK ( CHD! LEV 1 • EQ . 2 ) GO TO 76 

0040 


CBDILEV) - 2 
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FORTRAN* IV-P V02-51 19:03*28 04-APR-85 

CORSRT.FTN /TR:DLOCKS/VR 

0049 

0050 

005 1 76 

0052 

0053 

0054 


BDRY< 2 «LEV ) • DDRY<2,LEV> ♦ 1 

GO TO 72 

LEV - LEV - 1 

IF < LEV .GT . 0 ) GO TO 75 

RETURN 

END 
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FORTRAN~IV»P 

£HOFLI.FTN J 

0001 


V02-51 

/TR : BLOCKS/VR 


19)01:54 


04-APR-85 


1 


0002 

000 :? 

0004 

0005 

0006 


0007 

0000 

0009 

0010 
001 1 
0012 

0013 

0014 


0015 

0016 

0017 

0018 


0019 

0020 
002 1 
0022 

0023 

0024 

0025 

0026 
0027 
0020 

0029 

0030 

0031 

0032 

0033 

003 4 

0035 

0036 

0037 
0030 


C 

C 

C 


C 

C 

c 


SUQROUT INF ENOFL I < RPER . NP IXS , ISLS , ISSS ,EL,ES,NL,LINE, 

1 OLCNTR , NECNTR , MOLC , MEVC , MUM , OBEG I N , NBEG I M , 

2 OEND, MEND , ABUF ,PBUF ,L IM IT ) 


SEARCH FOR FINISHED OBJECTS 
IMTEGERM NPIXS< 1 > 

INTEGER OtCNTR< 1 >,HECNTR< 1 >,OBEGIN< 1 ),NBEGIN( 1 >,OEND< 1 ) ,NEND( 1 ) 
INTEGER I SLS< 1 > , I SSS< 1 ),EL< 1 ),ES< 1 ) 

REALM RPERi 1 >,ABUF< 1 >,PUUF( 1 ) 

IHTEGER-2 LIMITU, 1 ) 

SEARCH OLCNTR FOR I.D.’S MOT IN NECNTR 


C 

C 

c 


IF (L1NE.E0.HL) GO TO 52 
IF { MOLC . EG .0 ) RETURN 
IF < HE WC. EG. 0) GO TO 51 
1)0 10 J" 1 • NOI.C 
NEW ID - OLCNTR ( J ) 

DO 20 K-KNEVC 

IF < NE W I D.EU. NECNTR <tC ) ) GO TO 10 
20 CONTINUE 

DROPLET NOT CONTINUED 


C 

C 

C 


CALL PCHCAT< MOLC , OLCNTR, OBEG IN, OEND, NEWID, RPER > 

CALL F INAL < NEW ID,RPER,MPIXS, ISLS, ISSS , EL ,E5 , NUM ,ABUF , 
1 PBUF , L I M IT 1 

10 CONTINUE 

RETURN 

FINISH OFF EVERYTHING 

51 DO 50 J* 1 , NOLC 
MEW ID - OLCNTR ( 0 ) 

CALL PCNCAT < MOLC , OLCNTR ,OBEG IN , OEND , NEWID • RPFR > 

CALL F INAL < NEWID ,RPER, HP IXS, ISLS , ISSS , EL , E5 , NUM , ABUF , 
l POUF . LIMIT) 

50 CONTINUE 
RETURN 

52 IK < MOLC . EQ . 0 ) GO TO 60 
IP <NEWC.EO.tf> GO TO 51 
1)0 U/r J«* 1 , MOLC 

hewid - olcntrij) 

00 90 K»1,NEWC 

IF < NEWID. EQ . NECNTR! K ) ) CO TO 1 00 
90 CONTINUE 

CALL PCNCAT<NOLC, OLCNTR. OBEGIN, OEND. NEWID, RPER > 

CALL FINAL (NEWID. KPER.NP IXS, ISLS, ISSS. EL, E5, NUN, ABUF, 

1 POUF, LIMIT) 

GO TO 60 

100 CALL rCNCATlNLWC, NECNTR. NBEG I N , NEND , NEWID, RPER) 

CALL K INAL< NEWID. KPER.NP IXS. ISIS, ISSS. EL, LS. HUM, ABUF. 

1 PBUF, LIMIT) 

80 CONTINUE 

HE TURN 
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F0R7RAN fcV IV-P 

ENOFLI-FTN 


V02-51 

/7R : ULOCKS/WR 


19:01:54 


04-APR-85 


P 


0039 

0040 
004 1 

0042 

0043 

004 4 
004**i 
0046 


60 IF ( NEVC . EO . 0 ) RETURN 

do 70 j-i,mi:wc 

NEVID - MECNTR(J) 

CALL PCNCAT < NEVC , MECNTR , NREG IN ,NEND * NEVID , RPER ) 

CALL F I NAL ( NEVID , RPER , NP I XS , I SLS , I SSS , EL , ES , HUM , ABUF 
1 PBUF, LIMIT) 

70 CONTINUE 
RETURN 
END 



FORTRAN IV-PUJS V0 2-51 19:29:49 25-NAY-82 PAGF 1 

FlNAL.fTN /TR : BLOCKS/VR 


0001 SUBROUTINE F I NAL< NEWIO . RPER . NP I XS . I SLS ♦ ISSS . EL , ES , NUM .ABUF , 

1 PBUF ♦ L IM IT ) 

C 

C TO TERMINATE AN OBJECT 
C 

0002 IMTEGERM NPIXS(l) 

000:1 INTEGERS L IM IT< 4 , 1 ) 

0004 INTEGER I SL S< 1 ) , 1 3SS< 1) ♦ EL ( l ) , ES< I > 

0005 REALM RPER< 1 ),AUUF< 1 )*PBUF( 1 ) 

0H0U BYTE OR IGA< 512) 

0007 COMMON /MI MARE/ AMINAR 

0001) IF < NP I XS< NEW I B ) . EU .0 ) RETURN 

0000 MUM - NUM ♦ I 

0010 IF <NUM.LT. 201) GO TO 5 

D011 TYPE Number of 3pots 9 t buffer space available' 

0012 STOP 
C 

C ACCUMULATE DROPLET STATISTICS 
C 

0013 5 ABUF < MUM > • NPIXS(NEUID) 

0014 PBUF ( HUM ) - RPER< NEWII) ) 

0015 LIMIT! I, NUM) » ISI.S! MEVID I 

0010 L IM IT! 2 , HUM ) - ISSS<MEVID) 

0017 L (M IT< 3 « NUM ) - El.< NF.WID ) - ISLS(NEVID) ♦ 1 

00 1 0 L IM IT< 4 , MUM ) - ES<NEVID> - ISSS(NEUID) ♦ 1 

C CLEAR OUT BINS 

0010 IF < ABUF < HUM ).GT. AMINAR) GO TO 10 

0020 CALL DORRAR< 1, LIMIT! 1 , NUM ) ,L IMIT< 2, NUM ) ,L IMIT< 3, NUM ) «L IMIT< 4, NUM ) ) 

0021 HUM-NUM-1 
C 

0022 10 NP I X3< NEW ID ) - 0 

0023 RPER(HKWID) - 0. 

0024 ISI.S! NEVID) - 32000 

0025 1 SSS! NEWID ) - 32000 

0020 EM MEW ID) - 0 

0027 ES(HKWID) - 0 

002U RETURN 

0023 END 
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FORTRAN 'TV^P V02-51 10:50:30 Z3-FEB-84 P 1 

RAINST.FTN / TR :BLOCKS/VR 


0001 

C 

C 

C 

r; 

C 

0002 

D003 

0004 

0005 
0005 
000 / 

0000 

0000 

0010 

0011 

0012 

0014 

0014 

0015 
00 l 5 
00 17 
0014 
0014 
0020 
0021 

C 

C 

0022 451 

0024 

0024 

C 

0025 

0025 

0027 

0024 

0024 

0040 

004! 

0042 

0044 

0044 

0045 


0045 

0047 

4044 

0444 

0410 

0041 

0042 

0044 
004 4 

0045 


PROGRAM RA INST 

Road droplets* data flla and compute general statistic 

Programor: Miguel A. Neman 

REAL *4 DniJF< 1 200 ) , VOIJIJF ( 1200),ESFFRf 1200) 

REAL*!) SHU 4 f SMD2 , U'fOT , E37UT , SMD4 .COMUl) 

integer*2 home 

LOGICAL * 1 ANSA 

REAL *4 VOI.H1 S < 0i 74 >, FRENIS! 0:79 > 

REAL *4 AREA.PEIUM.D3.D2 J),U4 

LOG I CAL * 1 N020IH 12), NOZ I N < GO > , MOZ 10(00) 

E04 I VALENCE (H020IH 1 >,N02ID< 1 ) ) 

eou (valence? noz in? i i.nozid? i 3 > ) 

INTEGER *2 0PEPI.U 
LOG 1 CAL * 1 F ILL HA ( 22 > 

INTEGER *2 IllUFER? I) ) , 3L .S3, NL , N3 
EOU l VALENCE < IIHJFEU? l ),SL) 

EOO (VALENCIA IllUFER? 2 ) . SS ) 

EQUIVALENCE' MlUI EU? 4 > . HI. ) 

E OH I VALENCE? IllUFER? 4 1,115 ) 

E 0 4 I V A L E MC E (I R IJ F E I? ( 5 ) . All E A ) 

E Oil I VALENCE? IISUFER? 7 ) .PER IN 1 
DATA VOL 11 IS/UD*0. / , Fit EH 13/80*0./ 

DATA NQZUD / * H ' . ’O’. *2', 'Z'.'L'.'E', 1 I ',* D '/ 

Open data file 
FORMAT? 0.22A1 ) 

TYPE *,* Enter filename* 

ACCEPT 45 1 , LONGO , ? F I LENA? I ) , I - 1 ,LOHGO ) 

OPEN? UN IT-3 .HAHE-F 1LENA .TYPE* 'OLD ' . 

1 1) I 3P- * KEEP 1 .ACCESS- 'DIRECT 1 , RECORDS IZE-4 , FORM- ’ UNFORMATTED' ) 
READ? 3' 11 IIIUEER 
HUME- IllUFER? 1 > 

1COMT-0 

3HU4»0. 

SHI>4«0. 

3HD2-0. 

DTOT-0. 

ESTOT-IJ. 

PI-3. 14159 

QI-4 . /< 3 . *SORT< P I ) ) 

Main loop 
DO 100 K»2 , NOME 
READ? 4* K) IDUFER 
IE ? IllUFER? 1 >.CU,0l GO TO 200 
ICQNT-lCONT+l 
l)2-?4.*ARKA)/PI 
D-SORT? D2 ) 

113-1)2*1) 

1)4-112*1)2 

ESEEU? I COIIT >-? PER IM*PER IN )/< AREA*4 * *P I ) 

DDIJKI ICON T ) a D 


-Bl- 



10:50*30 


23-FEB-Q4 


FORtRAM 

"IV-P 

V02-51 10:50*30 23-FEB-Q4 

RAlNST. 

FTN 

/TIUBLOCKS/VR 

0040 


VnilUFi ICONT )®QI*!AREA*M*5) 

004/ 


3M»4-SNI)4-H>4 

004 il 


SIID3-3MI>3«M>3 

0040 


3HU2«3MI>2*D2 

0050 


IIT0T-DT0T4D 

0051 


E STUT **E STOT +E SF E R ! ICONT ) 

005.'! 


I I» 10*1) 

005 1 


l t »HI 110(79, 11 1 

0054 


Kin* il i si 1 1 i«fkciiis! n )*i 

0055 

1 00 

VOI.lt I S! I I >»VOLIIIS! I 1 >+U3 

0055 

CONTINUE 

0057 


CLOSE (UNIT-3) 

0O50 


uo to 300 

0050 

200 

CLOSE (UNIT-3 > 

0000 


TYPE Flic records counter error' 

006 1 

C 

0 

*: 

300 

CALL EXIT 


Begin general statistic 

0002 

TYPE *,• Enter no* 2 le code* 

0003 


ACCEPT 750, ITEX , ( N02 1N( 1 >, 1-1 , ITEX) 

0O04 

750 

FORMAT { 0 , 0I1A 1 ) 

0055 


TYPE Enter water pressure !psl)' 

0000 


ACCEPT * ,WATP 

000/ 


TYPE ' Enter wind tunnel speed ! ft/sec >' 

00011 


ACCEPT *,WTS 

0000 


TYPE Enter downstream coordinate (Inch)' 

0070 


ACCEPT * , DOWI. 

00/1 


PRINT 4 , H020U , ( NOZ IN! I >, I«1 , ITEX) 

00/2 


Pll I NT 5 ,WATP 

00/3 


PIIIMT li.WTS 

00/4 


PR1HT 7.DOWL 

00/5 


PIIIMT 11, ICOHT 

00/0 

4 

FORMAT < X ' 12A 1 , GOA 1 ) 

007/ 

5 

FORMAT! • Water pressure • , «F0*2'' psl*) 

0070 

0 

FORMAT!' Wind tunnel speed ••’♦F6.2, 1 ft/sec ’ ) 

0070 

7 

FORMAT! * Downstream location =*',FG.2'' Inch') 

0000 

11 

FORMAT! * Humber of Drops *',I5) 

0001 


PIIIMT 12 

0002 

12 

FORMAT! ' Length diameter <D10)') 

0003 


HAMO* ICONT 

0004 


CALL STATS! NANO , DBUF « DTOT ) 

0005 


PR HIT 13* SORT! SMD2/ FLOAT! ICONT ) ) 

0000 

13 

FORMAT!' Area moan diameter ! D20 )■ ’ « F 1 5 . 6 ) 

000/ 


C0H0I)**SND3*P I/O. 

0000 


PRINT 17 ,C0MUI1 

0000 

17 

FORMAT! ' Drops total vo 1 ume*» ' * F 15 . 6 ) 

0000 


HAMO* ICOHT 

000 1 


CALL STATS! HAMO , VQBUF , COMOD ) 

0002 


VOL ME l>“ SMI) 3 /FLOAT < ICOHT ) 

0003 


VULHED^VUI.HED**! 1 ./3. ) 

0M04 


PRINT 1 4 , VOL MED 

0005 

14 

FORMAT! * Volume mean diameter < D30)» ' « F 15 • 6 ) 

000O 


PRINT 10 « SMD3/3ND2 

000/ 

10 

FORMAT!' Sauter mean dlamoter ! D 32 )* ' , F 1 5 . 6 ) 

0000 


PRINT 30« SMD4/SHU3 
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FOUTRAII* IV^P V02-51 10:50:30 23-FEB-04 P 3 

RA INST ,FTN /TICULOCKS/WR 


0099 

30 

0100 


0101 

31 

0102 


0103 


0101 


0105 


0106 

008 

0107 


01011 


0109 


01 10 


0 111 


0 112 


0111 


0 114 

500 

0 115 


0 1 1 6 


0117 


01 1« 


0119 


0 1 20 


0121 


0122 


0 121 


0 12 1 

600 

0125 


0126 


0 1 2 / 


0120 


0129 


0130 


0131 

400 

0132 


0 133 


0 134 


0135 


0136 


0137 


013*1 


0139 


0140 

650 

014 1 

000 

0142 



FORMAT! ' Volume distribution mean diameter ( D43 ' , F 15 - 6 ) 

HR HIT 31 

FORMAT f ' Circularity (per I mate r * * 2 /area* 4 p 1 > ’ ) 

IIAHU-ICUNT 

CAM. STATS ( NANU , ESFKR , ESTOT ) 

TYPE * t * Not the d I str I but Ions <Y/HI' 

ACCEPT 0011, ANSA 
FORMAT < A 1 ) 

IF < AH5A.HE . ' Y * ) GO TO 000 
UPEPLO-0 

:;ahp«fi.oat< icomt)/100. 

SHD'!»SHI)J/10ir. 

DO 500 K-Jf./'J 
VOUIlSUd-VOLIIISUC I /SHOT 
FKEII I SI K >-l KEII1S! 1C l/SAMP 
CONTINUE 

TYPE Fro«i»ieiicy distribution <Y/M)' 

ACCEPT U0!t,AN5A 

IF < ANSA . HE . Tl 00 TO 600 

CALL RA INII I < FREII IS .OPEPLO ,0. , • i> 

IF < OPEPLO . EO.0 ) GO TO 000 
OPEPLO— OPEPLO 

CALL SYMUOU-4.3, l . 15, . 1, 'DROPS DIAMETER FREQUENCY DISTRIBUTION * « 

I 30., 37) 

CALL SYMBOL! - 1 . 35 , - . 75 , . 1 , * FREQUENCY < X > ' , 1 00. , 1 3 > 

00 TO 4110 

TYPE *,' Cumulative volume distribution < Y/N ) ' 

ACCEPT 1100 , ANSA 

IF (ANSA. HE. 'Y* ) 00 TO 650 

CALL RAIMIII! VOLIIIS. OPEPLO ,0. , . 1 ) 

IF ! OPEPLO. GE. 01 00 TO 650 

CALL SYMBOL! -4 . 3 , 1 .45 , . 1 , 'VOLUME CONTRIBUTION VS D IAMETER ’ , 90. ,31 } 
CALL SYMBOL!-. 05,-. 75. . 1 , 'VOLUME CONTRIBUTION < X ) 1 , 1 00 .,231 
CALL SYMBOL! .7,2. 35, . 1 , 'DIAMETER { MM ) ' , 90. , 1 3 ) 

CALL SYMBOL! -5.5,0. , . 1 ,110210,90. , IT E X + 1 2 ) 

CALL SYMBOL! -5.25,0. , . 1 , 'WATER PRESSURE ! PS I ) ' , 90. , 20 ) 

CALL NUMBER! -5.25 ,2.25, . 1 ,WATP ,90. ,-l ) 

CALL SYMBOL! -5 .0,0. , . 1 , ' WIND TUNNEL SPEED ! FT/SEC 90. , 26 ) 

CALL HUMBER! -5.0,2. 05, . 1,VTS,90. , -1 ) 

CALL SYMBOL! -4. 75, 0. I, 'DOWNSTREAM STATION ! INCH ) ' , 90. , 25 ) 

CALL HUMBER! -4.75.2.75, . 1 . DOVL ,90. ,-l ) 

II (OPEPLO. OT.0) 00 TO 600 

IF (OPEPLO. ME. 0) CALL PLOT! 0. ,0. , 999 ) 

CALL EXIT 
END 
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